NONPERTURBATIVE QCD WITH MODERN TOOLS 



C. D. ROBERTS 
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In these lectures I introduce and explore a range of topics of contemporary interest 
in hadronic physics: from what drives the formation of a nonzero quark condensate 
to the effect that mechanism has on light and heavy meson form factors, and 
the properties of the quark-gluon plasma. The trail leads naturally through a 
discussion of confinement, dynamical chiral symmetry breaking and bound state 
structure: phenomena that require nonperturbative methods for their explanation. 
In all of this, the necessary momentum-dependent modification of the quark and 
gluon propagators plays a significant role. 



1 Hadron Physics 

The strong interaction spectrum has a regularity that can be understood sim- 
ply. Mesons are composed of a constituent quark and its antiparticle, and 
baryons of three constituent quarks. The quarks have quantum numbers: 
colour - red, blue, green; flavour - up, down, strange, etc., and spin - up 
or down. Attributing a flavour-dependent mass to these constituent quarks: 
M u w M d ~ 320 MeV, M s ~ 400 MeV, and demanding that only colour sin- 
glets exist, one obtains a quantitative description that includes the multiplct 
structure. Almost. 

In this picture the pion is a pseudoscalar meson with the spins of its 
constituents antiparallel: ss 140 MeV. Align the spins and one has a vector 
meson, which is identified with the p: m„ ss 770 MeV. The simplest explanation 
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is hyperfine splitting, akin to the splitting of each level in the hydrogen atom 
that measures whether the spin of the electron and proton are aligned or not. 
A difficulty with this interpretation, however, is that the pseudovector- to 
scalar- meson mass ratio, w^/tn^, should have the same origin but it is much 
smaller. Equation (1) also shows that the excitation energy in the pseudoscalar 
channel, which is measured by the mass difference between the pion and its 
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first excited state, is much larger than that in the vector channel. As a final 
illustration of the peculiarity of the pion, I note that the expected value of the 
ratio in the third column is sa 2 because it should just count the number of 
constituents in a baryon as compared to a meson. It is for the p but not for the 
7r. The "unnaturally" low mass of the 7r has many consequences; e.g., it means 
that attraction in the N-N interaction persists over a much longer range than 
repulsion, which underlies the line of nuclear stability and makes possible the 
existence of heavy elements. 

Although the tt was the first meson discovered, its peculiar properties are 
still not understood widely. They are almost completely determined by the 
phenomena driving dynamical chiral symmetry breaking [DCSB]. That is an 
intrinsically nonperturbative effect, the understanding of which is one focus of 
these lecture notes. 

Quarks provide a means of understanding much of the regularity in the 
hadron spectrum. However, they do not form part of that spectrum: quarks 
carry the colour quantum number and the spectrum consists only of colourless 
objects. Are they then a mathematical artifice useful only as a means of realis- 
ing group theory? No, they are observed in inclusive, deep inelastic scattering: 
ep — ► e! + "debris" , even though the debris never contains an isolated quark. 

The cross section for deep inelastic scattering from a proton target is 

^ = iWI ( 2W ' lQ2 - "» 8i " 2 i + w ' {Q "- "> cos2 • (2) 

where 9 is the scattering angle, Q 2 is the spacelike squared momentum transfer, 
and v = E e > — E e for p p — 0. It is characterised by two scalar functions: Wi, 
W2, which contain a great deal of information about proton structure. If 
the proton is a composite particle composed of pointlike constituents, then it 
follows that for Q 2 — > 00 , v — > 00 with v Q 2 

Wi(Q 2 ,v 2 ) -> F 1 {Q 2 lv 2 ) , 

vW 2 (Q 2 ,v 2 )^F 2 (Q 2 /is 2 ); W 

i.e., that in inclusive, deep inelastic processes W1/2 are functions only of x := 
Q 2 jv 2 and not of Q 2 and v 2 independently. 

Such experiments were first performed at SLAC in the late sixties and. have 
been extended widely since then. Equation (||) is confirmed completelyp it is 
correct as Q 2 varies over four orders-of-magnitude. The presence of pointlike 
constituents is the only possible explanation of this behaviour, and many other 
observations made in such experiments. Hence quarks exist but are "confined" 
in colourless bound states, and understanding the origin and meaning of that 
confinement is another focus of these notes. 
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The role that quarks play in hadron structure is a primary subject of con- 
temporary experimental and theoretical nuclear physics and it is not necessary 
to employ high energy or large momentum transfer to explore it. The electro- 
magnetic form factor of the neutron is a good example. The differential cross 
section for electron-neutron scattering: 



da da 
dn = dQ 2 



Mott 



G 2 E (Q 2 ) + rG 2 M (Q 2 ) +2rG2M(Q2)tan2 _j (4) 



t := Q 2 /(4M^r), involves the "Mott" -term, which describes electron scattering 
from a point source, and the electric and magnetic form factors, Ge and Gm, 
that describe the structure of the neutron. Of particular interest is the electric 
form factor. For a neutral particle Ge{Q 2 = 0) = 0, and it is identically zero 
for a neutral point particle. Therefore deviations from zero are a good measure 
of ones description of the internal structure of thejieutron. 

The neutron charge radius has been measuredcl in the transmission of low 
energy neutrons through Pb atoms: 



a R dG E (Q 2 



= -0. 113 ± 0.003 fm 2 . (5) 

? 2 =o 



dQ 2 

It is negative and that must be understood. In studies&til of K° and K°* 
mesons one finds r 2 K0 < and r 2 K0 , < 0. These results are understood simply: 
the mesons are s-d bound states and, being lighter, the negatively charged 
d-quark is able to propagate further from the system's centre-of-mass thereby 
generating a negative charge-distribution at long-range. That mechanism gen- 
erates 60% of the observed value:B r\ = — 0.054 ± 0.026 fm 2 , and something 
similar is likely to be responsible for a large part of r 2 . 

The Q 2 - ev ahrtion f the form factors is more difficult to measure. In the 
absence of free neutron targets it must be inferred from the scattering of elec- 
trons from loosely-bound collections of polarised nucleons; e.g., deuterium or 
3 He. However, even in such loosely bound systems the inferred cross section 
is sensitive to details of the model used to calculate the nuclear wave function 
and the current operator, although it is anticipated that these effects can be 
minimised in double polarisation experimentstj The present status of measure- 
ments of G E is illustrated in Fig. |l|. 

Such complications often arise in prospecting for quarks and gluons. How- 
ever, their role in hadron structure is identifiable only in the signatures they 
write in, for example, the Q 2 -dependence of form factors, so the difficulties 
must be overcome. The accurate calculation of form factors on the entire 
domain of accessible Q 2 can be a very good test of any given description of 
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Figure 1: G^, - summary of dataQ. The lines are two-parameter fits using different TV- TV 
potentials in the calculation of the deuterium wave function. The "Q 2 -range" is that of a 
recent MAM1 [Mainz] experiment: 3 He(e, e'n). 



hadron structure. It relies on combining an understanding of confinement and 
DCSB with that of the intrinsically nonperturbative analysis of bound state 
structure. This synthesis is the primary focus of these notes. 

Given that at low-energy and small-Q 2 the strong interaction is charac- 
terised by the nonperturbative phenomena of confinement, DCSB and the be- 
haviour of colourless bound states, is its study ever simple? Yes, in scattering 
processes at high-energy or large-Q 2 where perturbation theory is applicable 



because QCD is asymptotically free. As described in Sec. 1.2, the coupling 
constant evolves with Q 2 : 

as(Q 2 ) Q! ^°°0, (6) 

and hence quarks and gluons behave as weakly interacting, massless particles 
when the momentum transfer is large. 

This also leads to the identification of another nonperturbative phenomenon. 
If one introduces the intensive variables: temperature, T, and quark chemical 
potential, /i, QCD acquires additional mass-scales with which the coupling can 
run. Hence, 



a S {Q 2 = 0; T,/x) ~ for T > A QCD and/or fj, > Aqcd 



(7) 
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where Aqcd ~ 200 MeV is the intrinsic, renormalisation-induced mass-scale in 
QCD. It follows that at finite-(T, /i) there is a phase of QCD in which quarks 
and gluons are weakly interacting, irrespective of the momentum transfer ;Q i.e., 
a quark-gluon plasma [QGP] . In this phase confinement and DCSB are absent 
and the nature of the strong interaction spectrum is qualitatively different. 
Phase transitions are intrinsically nonperturbative phenomena. 

1.1 Modern Experimental Tools 

Understanding the strong interaction, applying nonperturbative methods in 
QCD, is a difficult problem. As in all of physics, progress here and the devel- 
opment of understanding is promoted by the challenge of explaining real data, 
and there is a new generation of experimental facilities whose goal is to furnish 
this challenge. 

Deutsches Elektronen-Synchroton] in Hamburg is the site of the 
Hadron- Electron Ring Accelerator facility [HERA]. It is a collider: in a 6 336 m 
underground tunnel a longitudinally-polarised, 30 GeV er or e + beam collides 
with a counter-circulating, 820 GeV proton beam at four interaction zones. 
There are currently three large experimental collaborations using HERA: HI 
[since 1992], ZEUS [since 1992] and HERMES [since 1995], with the goal of 
elucidating the internal structure of the nucleon. 

Newport News, Virginia, is the site of the Thomas Jefferson National Ac- 
celerator Facility^ [TJNAF, formerly CEBAF]. This machine currently acceler- 
ates electrons to 4 GeV in five circuits of an oval-shaped, 1 400m underground 
tunnel. The electrons are then focused into three large, underground target 
areas: Halls - A [50m diameter], B [30m] and C [50m]. Although Hall B is not 
yet fully operational, experiments have been underway in Halls A and C for 
two years. Many have been completed and other important experiments are 
underway; e.g., Hall C is currently hosting a measurement of the electromag- 
netic pion form factor, which will probe the evolution from the nonperturbative 
to the perturbative domain in QCD, and that of is scheduled. The facility 
plans to expand the research capability of the accelerator by raising the beam 
energy to 12 GeV by the year 2006. 

Cornell University in Ithaca, New York, is the site of the Cornell Electron 
Storage Ringl3 [CESR]. A circular e + -e~ collider, 12m underground with a 
circumference of 768 m, CESR is a 6-quark factory. With the CESR collision 
energy tuned to that of the T(4s) 6-6-resonance, which is massive enough to 
decay to a B-B meson pair, the CLEO detector is able to observe the decays 
of -B-mesons. A careful analysis of these decays can provide tight constraints 
on the elements of the Cabibbo-Kobayashi-Maskawa [CKM] matrix, which 
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describes the mixing of the eigenstates of the strong and weak interactions. 
Most of what is currently known about i3-meson decays is due to experiments 
using the CLEO detector. However, since experiments using B-mesons are a 
far-reaching tool to probe the standard model, it is anticipated that the next 
five years will see the construction of three more 6-factories: at KEK in Japan; 
at SLAC; and in the completion of the HERA-B site at DESY. Heavy-quark 
physics is a burgeoning, contemporan£_.focus. 

Brookhaven National Laboratorytil [BNL] is the site of the AGS [Alternat- 
ing Gradient Synchrotron] and the future site of RHIC, the Relativistic Heavy 
Ion Collider. Its relativistic heavy ion programme currently focuses on the 
study of nuclear matter at extremes of temperature and density, and prospect- 
ing for the quark gluon plasma. The AGS has been operating since 1962 and 
its fixed target experiments are laying the foundations of the search for the 
QGP. In the future, the AGS will act as an injector for RHIC, which is due 
for completion in 1999. Using two concentric rings, each 3.8 km in diameter, 
RHIC will collide two 100 A GeV 197 Au beams at several detector sites around 
the ring. Producing a total centre-of-mass energy of ~ 40 TeV at each collision 
site, RHIC is expected to produce an equilibrated QGP, approaching it along 
a high-T, baryon-poor trajectory. 

Also, in the hunt for the QGP is the SpS [Super-proton Synchrotron] at 
CERN[L3 seven kilometres in circumference. Fixed target experiments with 
S+Pb at ~ 200 A GeV and Pb + Pb at ~ 160 A GeV are possible using this 
machine, producing much higher centre-of-mass energies than those available 
at the AGS. With these resources, experiments CERES/NA45 and NA50 have 
observed interesting effects: an enhancement in the dilepton spectrum and 
a suppression of the '(/'-production cross section. Whether these signal the 
formation of a QGP is an unresolved issue but their explanation is certainly 
providing a challenge to theory. The next phase is the construction of the LHC 
[Large Hadron Collider]. Currently anticipated to be completed in 2005, it 
will share the 27 km LEP tunnel and be capable of hosting heavy ion collisions 
that generate a centre-of-mass energy thirty-times greater than that possible 
at RHIC. If the QGP is not discovered with RHIC then the LHC can still 
promise to do so. If it is, then the LHC should provide a tool for thoroughly 
exploring its properties, and for many other things. 

I have surveyed only a few of the many existing and planned experimental 
facilities. Nevertheless, the intensity of effort being expended in prospecting 
for quarks and gluons, and in elucidating the subtle properties of QCD, is 
obvious in this. More detailed information is available at the web sites I have 
indicated and obvious searches will yield details about other facilities. 
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1.2 Perturbative QCD 

QCD is a gauge theory, which in some respects is like QED. The elementary 
fermions, quarks in this case, interact by exchanging massless vector bosons: 
gluons. In QCD the coupling is through the colour quantum number: the 
colour charge is the analogue of the electric charge. The special feature of 
QCD is that the gluon carries the colour charge, as well as the quarks, and that 
means the gluons are self-interacting, which is enough to completely change 
the qualitative character of the theory. 

In both QED and QCD the coupling "constant" is not constant at all but 
evolves with Q 2 and it can be calculated perturbativcly: 

a 1 . . 

1 — a/37rln {Q^/m^) 137 

= 7 r, (9) 

(HAT c -2iV / )ln(Q 2 /A2 CD J 

where the "lliV c ", N c = 3 is the number of colours, appears because the glu- 
ons are self-interacting. Nf is the number of fermion flavours. As depicted in 
Fig. ||, as(Q) decreases with increasing-Q, which is opposite to the behaviour 
of Q!qed(Q) quoted in Eq. (||). The rate of the evolution is also very differ- 
ent, being determined by vastly different mass-scales: Aqcd/toe ~ 450. The 
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reduction in as with increasing Q is called "asymptotic freedom" . It is due to 
gluon self-interactions and entails that perturbation theory is valid at large-Q 2 
in QCD; i.e., at short-distances, not at long-distances as in QED. 

A very successful application of perturbative QCD is the calculation of 

j£ ^e+e — — >hadrons (10) 

which is known to order a^: 



\ " 2 



R = N c ^ e f 
/=i 



1 + + 1.411 f^) 2 - 12.8 f^) : 



(11) 



where e/ is the electric charge of a quark of flavour /.|^| The theoretical pre- 
diction is compared with experiment in Fig. |^. It proves that N c = 3 because 
the normalisation of the prediction is otherwise incorrect. The figure also 
highlights the current-quark mass thresholds and confirms the electric charge 
assignments: below the c-quark threshold R s» iV c J2jf^ 3 = 2, then R rj 

N c X^/^i -4 e / = 10/3, an d above the 6-quark threshold R w N c Y^Li^ e "] = 
11/3. 

Perturbation theory is the most widely used, systematic tool in physics 
and in QCD it is extremely accurate in high-energy processes. In addition to 
the calculation of R, there are many other applications which confirm QCD 
as the theory describing the strong interaction. However, there are domains 
and problems in QCD that perturbation theory simply cannot describe. This 
is very obvious in Fig. || where the coupling is seen to increase with the inter- 
particle separation: QCD becomes a strong coupling theory for Q < 2 GeV 
[x > 0.1 fm]. Further, as observed in Sec. |l|, confinement, DCSB, bound state 
structure and phase transitions are intrinsically nonperturbative phenomena. 
They are at the core of hadron physics and hence nonperturbative methods 
are essential in prospecting for quarks and gluons. 



1.3 Nonperturbative Effects in QCD 

The extraordinary effects in QCD can be explained in terms of the proper- 
ties of dress ed-quark and -gluon propagators. They describe the "in-medium" 
propagation characteristics of QCD's elementary quanta, with the "medium" 
being the nontrivial ground state of QCD. A photon propagating through a 

"This expression assumes that all the quarks are massless but it is also known for massive 
quarks. In comparing with data, b-quark mass corrections are important. 
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Figure 3: Selected measurements of RU The positions of the ip, iji{2S) and four lowest 
T(nS) resonances are also indicated, which are not predicted by perturbation theory. A fit 
to the data in the lower panel yields a<j(Q = 34 GeV) = 0.142 ± 0.03. The experimental 
error is larger than the theoretical uncertainty. The plotted calculations include electroweak 
corrections, which predict the increase a,s liicm -♦ M z . 



dense e~-gas provides a familiar example of the effect a medium has on the 
propagation of an elementary particle. Due to particle-hole excitations the 
propagation of the photon is modified: 

Q 2 + ml' (12) 

i.e., the photon acquires an effective mass. The "Debye" mass, oc hp, 
the Fermi momentum, and it screens the interaction so that in the dense e~- 
gas the Coulomb interaction has a finite range: r oc 1/mu- Quark and gluon 
propagators are modified in a similar way. They acquire momentum- dependent 
effective masses, which have observable effects on hadron properties. 

1.3.1 Chiral symmetry. Gauge theories with massless fermions have a chiral 
symmetry. Its effect can be visualised by considering the helicity: A oc J -p, the 
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projection of the fcrmion's spin onto its direction of motion. A is a Poincarc 
invariant spin observable that takes a value of ±1. The chirality operator can 
be realised as a 4 x 4-matrix, 75, and a chiral transformation is then represented 
as a rotation of the 4 x 1-matrix quark spinor field 

q(x) ^ e^ e q(x) . (13) 

A chiral rotation through 9 = ir/2 has no effect on a A = +1 quark, q\=+ — ► 
q\— + , but changes the sign of a A = —1 quark field, q\—- — > — q\=-- In 
composite particles this is manifest as a flip in their parity: J p=+ <-> J p= ~; 
i.e., a 8 = tt/2 chiral rotation is equivalent to a parity transformation. Exact 
chiral symmetry therefore entails that degenerate parity multiplets must be 
present in the spectrum of the theory. 

For many reasons, the masses of the u- and d-quarks are expected to be 
very small; i.e., m u ~ ma <§; Aqcd- Therefore chiral symmetry should only be 
weakly broken, with the strong interaction spectrum exhibiting nearly degen- 



erate parity partners. The experimental comparison is presented in Eq. (14): 

770) 



iV(i + ,938) 



' 2 _ 
N(i ,1535) 



tt(0-,140) 
a (0+,980) 



(14) 
oi(l+,1260) 



Clearly the expectation is very badly violated, with the splitting much too 
large to be described by the small current-quark masses. What is wrong? 

Chiral symmetry can be related to properties of the quark propagator, 
S(p). For a free quark 

m-i-y-p 

So{p) = — 2 : 2 , 15 
m z + p z 

with {7^, /i — 1, . . . , 4} the Dirac matrices, and as a matrix 

S (p) - e^ s S (p)e^ e = + e^ e (16) 

p z + p z +m z 

under a chiral transformation. As anticipated, for m = 0, So(p) — > So(p); i.e., 
the symmetry breaking term is proportional to the current-quark mass and it 
can be measured by the "quark condensate" 

/d^p f d^p vn 

j^tv[S(p)}ocJ —- r — 2 , (17) 

which is the "Cooper-pair" density in QCD. For a free quark the condensate 
vanishes if m = but what is the effect of interactions? 
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Figure 4: A rotationally invariant but unstable extremum of the Hamiltonian obtained with 
the potential V(a, vr) = (cr 2 + tt 2 - l) 2 . 



Interactions dress the quark propagator so that it takes the form 

Sh) •= 1 = -H-PA(P 2 ) + B( P 2) 
W ' z 7 -p+£(p) p 2 4 2 (p 2 ) + B 2 (p 2 ) ' { ' 

where £(p) is the self energy, expressed in terms of the scalar functions: A and 
B, which are p 2 -dependent because the interaction is momentum-dependent. 
On the valid [weak-coupling] domain they can be calculated in perturbation 
theory and at one-loop order 



B(p 2 ) = m[ 1-^ln 

4-7T 



P 2 



(19) 



which is oc to. This result persists: at every order in perturbation theory every 
mass-like correction to S(p) is oc to so that m is apparently the only source of 
chiral symmetry breaking and (qq) oc m — > as m - > 0. The current-quark 
masses are the only explicit chiral symmetry breaking terms in QCD. 

However, symmetries can be "dynamically" broken. Consider a point- 
particle in a rotationally invariant potential V(cr, tt) = (a 1 + tt 2 — l) 2 , where 
(cr, tt) are the particle's coordinates. In the state depicted in Fig. [|, the particle 
is stationary at an extremum of the action that is rotationally invariant but 
unstable. In the ground state of the system, the particle is stationary at any 
point (cr, tt) in the trough of the potential, for which a 2 + tt 2 — 1. There are 
an uncountable infinity of such vacua, \9), which are related one to another by 
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Figure 5: DSE for the dressed-quark self-energy. The kernel of this equation is constructed 
from the dressed-gluon propagator (D - spring) and the dressed-quark-gluon vertex (r - open 
circle). One of the vertices is bare (labelled by 7) as required to avoid over-counting. 

rotations in the (a, 7r)-plane. The vacua are degenerate but not rotationally in- 
variant and hence, in general, (0\cr\6) =^ 0. In this case the rotational invariance 
of the Hamiltonian is not exhibited in any single ground state: the symmetry 
is dynamically broken with interactions being responsible for (8\a\9) =/= 0. 

1.3.2 Dynamical chiral symmetry breaking. The analogue in QCD is 
(qq) when m = 0. At any finite order in perturbation theory that is 
impossible. However, using the Dyson- Schwingcr equation [DSE] for the quark 
self energy [the QCD "gap equation"]: 

ij-pA(p 2 ) +B(p 2 ) = Z 2 i 1 -p +Z 4 m (20) 

+Zl J (tyi 92 D ^ ip ~ i] 7 -T * 7 • IMF) + B{P) va ^ p) ' 

depicted in Fig. |5|, it is possible to sum infinitely many contributions]^ That 
allows one to expose effects in QCD which are inaccessible in perturbation 
theory. 

The quark DSE is a nonlinear integral equation for A and B and its non- 
linearity is what makes it possible to generate nonperturbative effects. The 
kernel of the equation is composed of the dressed-gluon propagator: 

^W-^-^ gW.g^^p^, (21) 

b In Eq. (pc|), m is the A-dependent current-quark bare mass and J represents mnemonically 
a translationally-invariant regularisation of the integral, with A the regularisation mass-scale. 
The final stage of any calculation is to remove the regularisation by taking the limit A — » 00. 
The quark-gluon- vertex and quark wave function renormalisation constants, Zi(£ 2 , A 2 ) and 
•^2(C 2 ;^ 2 )i depend on the renormalisation point, and the regularisation mass-scale, as 
does the mass renormalisation constant Z m (( 2 , A 2 ) := Z2(( 2 , A 2 ) _1 Z4((^ 2 , A 2 ). 
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Figure 6: Illustrative forms of G(Q) that agree with the perturbative result at Q > 1 GeV. 
Three possibilities are canvassed: Q(Q = 0) < 1; S{Q = 0) = 1; and Q(Q = 0) > 1. 

where II(fc 2 ) is the vacuum polarisation, which contains all the dynamical infor- 



One chooses an initial seed for the solution functions: qA and oB, and evaluates 
the integral on the right-hand-side (r.h.s.). The bare propagator values: $A = 1 
and qB — m are often adequate. This first iteration yields new functions: \A 
and i-B, which are reintroduced on the r.h.s. to yield 2A and 2B, etc. The 
procedure is repeated until n A = n +\A and n B = n +±B to the desired accuracy. 

It is now simple to illustrate DCSB. Using the bare vertex and G(Q) de- 
picted in Fig. ||, I solved the quark DSE in the chiral limit. If Q(Q = 0) < 1 
then B(p 2 ) = is the only solution. However, when Q(Q = 0) > 1 the equa- 
tion admits an energetically favoured B(p 2 ) ^ solution; i.e., if the coupling 
is large enough then even in the absence of a current-quark mass, contrary to 
Eq. (p^), the quark acquires a mass dynamically and hence 




(qq) oc 



/ 



d 4 p B(p 2 ) 



for m = . 



(23) 



(2tt) 4 p 2 A(p 2 ) 2 + B(p 2 ) 2 
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This identifies a mechanism for DCSB in quantum field theory. The 
nonzero condensate provides a new, dynamically generated mass-scale and if 
its magnitude is large enough [] it can explain the mass splitting between parity 
partners, and many other surprising phenomena in QCD. The simple model 
illustrates that DCSB is linked to the long-range behaviour of the fermion- 
fermion interaction and the same is true of confinement. The question is then: 
How does D^ v (k) behave in QCD? 



2 Dyson-Schwinger Equations 

In the last section I introduced the DSE for the quark self energy. It is one 
of an infinite tower of coupled integral equations, with the equation for a 
particular n-point function involving at least one m > n-point function; e.g., 
the quark DSE involved the dressed-gluon propagator, a 2-point function, and 
the dressed-quark-gluon vertex, a 3-point function. The collection of DSEs 
provide a Poincare invariant, continuum approach to solving quantum field 
theories. However, as an infinite collection of coupled equations, a tractable 
problem is only obtained if one truncates the system. Historically this has 
provided an impediment to the application of DSEs: a priori it can be difficult 
to judge whether a particular truncation scheme will yield qualitatively or 
quantitatively reliable results for the quantity sought. As we saw, the DSEs 
are integral equations and hence the analysis of observables is a numerical 
problem. Therefore a critical evaluation of truncation schemes often requires 
access to high-speed computers]^ With such tools now commonplace, this 
evaluation can be pursued fruitfully. 

The development of efficacious truncation schemes is not a purely numer- 
ical task, and neither is it always obviously systematic. For some, this last 
point diminishes the appeal of the approach. However, with growing commu- 
nity involvement and interest, the qualitatively robust results and intuitive 
understanding that the DSEs can provide is becoming clear. Indeed, those 
familiar with the application of DSEs in the late- 70s and earkfa&Cte might be 
surprised with the progress that has been made. It is now cleaniiiEJ that trun- 
cations which preserve the global symmetries of a theory; e.g., chiral symmetry 
in QCD, are relatively easy to define and implement and, while it is more dif- 
ficult to preserve-local gauge symmetries, much progress has been made with 
Abelian theories^ and more is being learnt about non-Abelian ones. 

The simplest truncation scheme for the DSEs is the weak-coupling expan- 

c — (qq) 1 ^ 3 need only be one order-of-magnitude larger than in u ~ mj. 

d The human and computational resources required are still modest compared with those 
consumed in contemporary numerical simulations of lattice-QCD. 
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Figure 7: DSE for the gluon vacuum polarisation and propagator: solid line - quark; spring - 
gluon; dotted-line - ghost. The open circles are irreducible vertices. As indicated, the quark 
loop acts to screen the charge, as in QED, while the gluon loop opposes this, anti-screening 
the charge and enhancing the interaction. 

sion, which shows that they contain perturbation theory; i.e, for any given 
theory the weak-coupling expansion generates all the diagrams of perturba- 
tion theory. However, as with DCSB, the most important feature of the DSEs 
is the antithesis of this weak-coupling expansion: the DSEs are intrinsically 
nonperturbative and their solution contains information that is not present in 
perturbation theory. They are ideal for studying the phenomena I identified 
as the core of these lectures and in this application they provide a means of 
elucidating identifiable signatures of the quark-gluon substructure of hadrons. 

2.1 Gluon Propagator 

In Landau gauge the dressed-gluon propagator has the form in Eq. (^) and 
satisfies the DSE [a nonlinear integral equation] depicted in Fig. [?]. A weak 
coupling expansion of the equation reproduces perturbation theory and shows 
directly that in the one-loop expression for the running coupling constant, 
Eq. the "lliV c " comes from the charge-antiscreening gluon loop and the 
u 2Nf" from the charge-screening fermion loop. This illustrates how the non- 
Abelian structure of QCD is responsible for asymptotic freedom and suggests 
that confinement is related to the importance of gluon self-interactions. 

Studies of the gluon DSE have been reported by many authorsta with the 
conclusion that if the ghost-loop is unimportant, then the charge-antiscreening 
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Figure 8: Dressed-quark mass-function obtained in solving the quark DSE. 



3-gluon vertex dominates and, relative to the free gauge boson propagator, the 
dressed gluon propagator is significantly enhanced in the vicinity of k 2 = 0. 
The enhancement persists to k 2 ~ 1-2 GeV 2 , where a perturbative analysis 
becomes quantitatively reliable. In the neighbourhood of k 2 = the enhance- 
ment can be representecO as a regularisation of 1/fc 4 as a distribution. A 
dressed-gluon propagator of this type generates confinement and DCSB with- 
out fine-tuning, as I will elucidate. 



2.2 Quark Propagator 

In a covariant gauge the dressed-quark propagator can be written in a number 
of equivalent forms: those in Eq. ([l8]) and also as 

S(p) = -ir, . P a v (p 2 ) + a s ( P 2 ) = ^ ^2 {p2) ■ (24) 

As depicted in Fig. |L solving the quark DSE, Eq. (^Q), using a dressed-gluon 
propagator of the type described above, one obtains a quark mass-function, 
M(p 2 ), that mirrors the enhancement of the dressed-gluon propagator. The 
results in the figure were obtainedll3 with current-quark masses corresponding 



1G 



to 



m lGcV m 1GeV m 1GcV m 1 GcV 
m u/d m s c b 

6.6 MeV 140 MeV 1.0 GeV 3.4 GeV. 1 ; 



The quark DSE was also solved in the chiral limit, which in QCD is ob- 
tained by setting the Lagrangian current-quark bare mass to zcroll2 One ob- 
serves immediately that the mass-function is nonzero even in this case. That 
is DCSB: a momentum-dependent quark mass, generated dynamically, in the 
absence of any term in the action that breaks chiral symmetry explicitly. This 
entails a nonzero value for the quark condensate in the chiral limit. That 
M(p 2 ) 7^ in the chiral limit is independent of the details of the infrared 
enhancement in the dressed-gluon propagator. 

Figure || illustrates that for light quarks (u, d and s) there are two distinct 
domains: perturbative and nonperturbative. In the perturbative domain the 
magnitude of M(p 2 ) is governed by the the current-quark mass. For p 2 < 
1 GeV 2 the mass-function rises sharply. This is the nonperturbative domain 
where the magnitude of M (p 2 ) is determined by the DCSB mechanism; i.e., the 
enhancement in the dressed-gluon propagator. This emphasises that DCSB is 
more than just a nonzero value of the quark condensate in the chiral limit! 

The solution of p 2 = M 2 {p 2 ) defines a Euclidean constituent-quark mass, 
M E . For a given quark flavour, the ratio £/ :— Mj? /rrvj is a single, quanti- 
tative measure of the importance of the DCSB mechanism in modifying the 
quark's propagation characteristics. As illustrated in Eq. (prj), 

(26) 



flavour 


u/d 


s 


c 


b 


t 


mC ~2() UcV 


150 


10 


2.3 


1.4 


-> 1 



this ratio provides for a natural classification of quarks as either light or heavy. 
For light-quarks £/ is characteristically 10-100 while for heavy-quarks it is only 
1-2. The values of £/ signal the existence of a characteristic DCSB mass-scale: 
M x . At p 2 > the propagation characteristics of a flavour with < M x are 

altered significantly by the DCSB mechanism, while for flavours with m'j 
M x it is irrelevant, and explicit chiral symmetry breaking dominates. It is 
apparent from Eq. (|26| ) that M x ~ 0.2 GeV ~ Aqcd- This foims a basis for 
many simplifications in the study of heavy-meson observablestj 



2.3 Confinement 

Confinement is the absence of quark and gluon production thresholds in colour- 
singlet-to-singlet 5-matrix amplitudes. That is ensured if the dressed-quark 
and -gluon propagators do not have a Lehmann representation. 
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To illustrate a Lehmann representation, consider the 2-point free-scalar 
propagator: A(fc 2 ) = l/[k 2 + m 2 ]. One can write 




oo 



da 



■ 

z + a 



(27) 



where in this case the spectral density is 



p( x ) '■— 7; — " nm IA(— x ~ * e ) — A(— £ + ie)\ = S(m 2 — x) , 



(28) 



which is non-negative. That is a Lehmann representation: each scalar function 
necessary to specify the n-point function completely has a spectral decompo- 
sition with non-negative spectral densities. Only those functions whose poles 
or branch points lie at timelike, real-fc 2 have a Lehmann representation. 

The existence of a Lehmann representation for a dressed-particle propaga- 
tor is necessary if the construction of asymptotic "in" and "out" states for the 
associated quanta is to proceed; i.e., it is necessary if these quanta are to prop- 
agate to a "detector" . In its absence there are no asymptotic states with the 
quantum numbers of the field whose propagation characteristics are described 
by the propagator. Structurally, the nonexistence of a Lehmann representation 
for the dressed-propagators of elementary fields ensures the absence of pinch 
singularities in loops and hence the absence of quark and gluon production 
thresholds. 

The mechanism can be generalised and applied to coloured bound-States, 
such as colour-antitriplet quark-quark composites (diquarks) . A studyO of the 
quark-quark scattering matrix shows that it does not have a spectral decom- 
position with non-negative spectral densities and hence there are no diquark 
bound states. The same argument that demonstrates the absence of diquarks in 
the spectrum of SU(N C = 3) also provesEJ that in SU(N C = 2) the "baryons" , 
which are necessarily diquarks in this theory, are degenerate with the mesons. 

Dressed-gluon propagators with the infrared enhancement described above 
do not have a Lehmann representation and using forms like this in the ker- 
nel of the quark DSE yields a dressed-quark propagator that also does not 
have a Lehmann representation. In this sense confinement breeds confinement, 
without fine-tuning. 

2.4 Hadrons: Bound States 

The properties of hadrons can be understood in terms of their substructure by 
studying covariant bound state equations: the Bethe-Salpeter equation [BSE] 
for mesons and the covariant Fadde'ev equation for baryons. The mesons have 
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Figure 9: First two orders in a systematic expansion of the quark-antiquark scattering kernel. 
In this expansion, the propagators arc dressed but the vertices are bare. 

been studied most extensively and their internal structure is described by a 
Bethe-Salpeter amplitude obtained as a solution of 



where X H{q;P) := S{q + )T H {q;P)S{q.); 5(g) = diag(S„(g), S d (q), S.(q), . . 



g_l_ = q + rjp P, g_ = g — (1 — rjp) P, with P the total momentum of the bound 
state; and r,. . . ,u represent colour-, Dirac- and flavour-matrix indices. The 
amplitude for a pseudoscalar bound state has the form 



where T is a flavour matrix that determines the channel under consideration; 



e.g., T K+ :=_(l/2) (A 4 + i\ 5 ), with {X> , j = 1 ... 8} the Gell-Mann matrices. 



In Eq. tegj), K is the renormalised, fully-amputated, quark-antiquark scat- 
tering kernel and important in the successful application of DSEs is that it has 
a systematic skeleton expansion in terms of the elementary, dressed-particle 
Schwinger functions; e.g., the dressed-quark and -gluon propagators. The ex- 
pansion introduced in Ref. |l| provides a means of constructing a kernel that, 
order-by-order in the number of vertices, ensures the preservation of vector 
and axial- vector Ward-Takahashi identities; i.e., current conservation. 

In any study of meson properties, one chooses a truncation for K. The 
BSE is then fully specified and straightforward to solve, yielding the bound 




(29) 



T H (k;P) = T H 75 iE H {k;P)+<y-PF H (k;P) 



(30) 



+ 7 -kk-PG H (k;P) + a IMU k^P v H H (k; P) 
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state mass and amplitude. The "ladder" truncation of K combined with the 
"rainbow" truncation of the quark DSE [T^ — * 7^ in Eq. (pC|)] is the simplest 
and most often used. The expansion of Fig. ^| provides the explanation!!] for 
why this Ward-Takahashi identity preserving truncation is accurate for flavour- 
nonsinglet pseudoscalar and vector mesons: there are cancellations between the 
higher-order diagrams. It also shows why it provides a poor approximation in 
the study of scalar mesons, where the higher-order terms do not cancel, and 
for flavour-singlet mesons, where it omits timclikc gluon exchange diagrams. 



3 A Bound State Mass Formula 

The dressed-axial-vector vertex satisfies a DSE whose kernel ia,K, and because 



of the systematic expansion described in Sec. 2.4 it followsES that the axial- 
vector Ward-Takahashi identity [AV-WTI]: 

-iPJ%(lr,P)= (31) 
5- 1 (fc+) 7 5^- + 76^-«S- 1 (*:-) " M (0 rf (fc; P) - rf (k; P) M (f) , 

[M/q = diag(m£, m^, , . . .) is the current-quark mass matrix] is satisfied in 
any thoughtful truncation of the DSEs. That entails many important results. 

fn- The axial- vector vertex has a pole at P 2 — —m 2 H whose residue is 
fa-, the lcptonic decay constant: 

f H P, = Z ^f-^l i tr [CO* l^S{q + )T H {q; P)S(q-)} , (32) 

with the trace over colour, Dirac and flavour indices. This expression is exact: 
the dependence of Zi on the renormalisation point, regularisation mass-scale 
and gauge parameter is just that necessary to ensure that the l.h.s. is inde- 
pendent of all these things. 

Goldstone's Theorem: In the chiral limit 

f H E H (k; 0) = B Q (k 2 ) , F R (k; 0) + 2 f H F H (k; 0) = A (k 2 ) , 

G R (k- 0) + 2 f H G H {k- 0) = 2A' (k 2 ) , H R (k- 0) + 2 f H H H (k; 0) = , 

(33) 

where A (k 2 ) and B (k 2 ) are the solutions of Eq. ( p0|) in the chiral limit, and 
Fri G r and H R are calculable functions in This shows that when chiral 
symmetry is dynamically broken: 1) the flavour- nonsinglet, pseudoscalar BSE 
has a massless solution; 2) the Bethe-Salpeter amplitude for the massless bound 
state has a term proportional to 75 alone, with the momentum-dependence 
of En(k;0) completely determined by that of Bo(k 2 ), in addition to terms 
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proportional to other pseudoscalar Dirac structures that are nonzero in general; 
and 3) the axial- vector vertex, T^(k;P), is dominated by the pseudoscalar 
bound state pole for P 2 ~ 0. The converse is also true. Hence, in the chiral 
limit, the pion is a massless composite of a quark and an antiquark, each 
of which has an effective mass M E ~ 450 MeV. Wit h an infrared enhanced 
dressed-gluon propagator of the type described in Sec. |2.l| , this occurs without 
fine-tuning. 

In-meson Condensate: The pseudoscalar vertex also has a pole at P 2 = 
—m 2 H whose residue is 

irn = Z,f^ ~tr [(T ff ) 4 ^S(q + )T H (q; P)S(q-)] := - ^ . (34) 

The renormalisation constant Z4 on the r.h.s. depends on the gauge parameter, 
the regularisation mass-scale and the renormalisation point. This dependence 
is exactly that required to ensure rn is finite in the limit A — > 00 and gauge- 
parameter independent, (qq)^ is the in-meson quark condensate. 

Mass Formula: There is an identity between the residues of the pseu- 
doscalar meson pole in the axial-vector and pseudoscalar vertices that is sat- 
isfied independent of the magnitude of the current-quark mass: 

f H m H = - M H (qq)f , M H := tr flavour [m (c) [t h , (T ff )'}] , (35) 

e.g., for the 7r: A4h = Tn^ + This is a mass formula for flavour-octet pseu- 
doscalar mesons and the r.h.s. does not involve a difference of massive quark 
propagators: a phenomenological assumption often employed. The renormali- 
sation point dependence of (qq)^ is exactly such that the r.h.s. of Eq. (|3^) is 
renormalisation point independent. 

3.1 A Corollary 

For small current-quark masses, using Eqs. d3(]) and (|33|), Eq. (S3) yields 



1 ' A 



r H = - W (??>C . " (5?>C := Z 4(C 2 , A 2 ) N c / tr Dirac [S 7% =o(q)} , (36) 

JH 



'I 

where the superscript "0" denotes that the quantity is evaluated in the chiral 
limit and (qq)®, as defined here, is the chiral limit vacuum quark condensate; 
i.e., the in-meson condensate becomes the vacuum quark condensate in the 
chiral limit. One obtains immediately from Eqs. (]35|) and ( ^6|) 

Uml = -\mi + mft (qq) , + O {m 2 ) , (37) 



d 

m. 



gj 

ik '".k - l' a i + mC s] (<7<7}° + (rh 2 q ) , (38) 
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which exemplify what is commonly known as the Gell-Mann-Oakes-Renner 
relation. \ rh„ is the renormalisation-point-independent current-quark mass.] 
Since Eq. (|35| ) is valid independent of the current-quark mass, it is also appli- 
cable to heavy pseudoscalar mesons. A corollary valid in this sector is described 
in Sec. ||. 

4 An Exemplar 

I have already made use of a modclEl to illustrate some of the robust results 
of DSE studies and I will now describe and explore that model in more detail. 
For the kernel of the quark DSE it employs the analogue of the lowest-order 
BSE kernel in Fig. |: 

Zi J ^3 2 ^(p-g)y7/,%)r^,p) (39) 

a 1 nll„ „\2\ nfree/ _\ A . o/„\ A . 



G((p - qY) D*™(p - g)_ 7p 5(g)— 7 , • 



(2tt) 4 ^ HJ ' ^ ^ H; 2 ,il yH ' 2 

This is the "rainbow" approximation, in which the specification of the model 
is complete once a form is chosen for the "effective coupling" Q(k 2 ), 

The requirement of consistency with the AV-WTI motivates the Ansatz: 

G{Q 2 ):=4ira{Q 2 ), (40) 

so that the form of G(Q 2 ) at large-Q 2 is fixed by that of the running coupling 
constant. This Ansatz is often described as the "Abelian approximation" be- 
cause the l.h.s. equals the r.h.s. in QED. In QCD, equality between the two 
sides of Eq. (EH) cannot be obtained easily by a selective resummation of dia- 
grams. As reviewed in Ref. flqj , Eqs. (5.1) to (5.8), it can only be achieved by 
enforcing equality between the renormalisation constants for the ghost-gluon 
vertex and ghost wave function: Z± = Z3. 

The explicit form of the Ansatz employed in Rcf. is depicted in Fig 
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and its qualitative features are easily understood. In the infrared it has an 
integrable singularityEHI oc S 4 (k) and a finite- width approximation to <5 4 (fc), 
normalised such that it has the same J d 4 k as the first term, and in the ul- 
traviolet it is dominated by a(k 2 )/k 2 . However, it does not have a singularity 
on the real-A: 2 axis, which ensures gluon confinement through the absence of a 



Lehmann representation, Sec. 2.3 



The model has ostensibly three parameters: D, a mass-scale, and u> and 
TO*, two range parameters. However, in the numerical studies the values lu — 
0.3GeV[= l/(-66fm)] and m t = 0.5GeV[= l/(.39fm)] were fixed, and only D 
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and the rcnormalised u/d- and s-current-quark masses varied in order to satisfy 
the goal of a good description of low-energy tt- and iC-meson properties. This 
was achieved with 



D = 0.781 GcV 2 , 



u/a 



3.74 MeV, m$ = 82.5 MeV 



(41) 



at £ w 20 GeV, which is large enough to be in the perturbative domain. 



As remarked in Sec. L3, the chiral limit is unambiguously defined by m = 
so there is no perturbative contribution to the scalar piece of the quark 
self energy, B(p 2 ); and in fact there is no scalar, mass-like divergence in the 
perturbative calculation of the self energy. 

Figure [ll] depicts the dressed-quark mass function, M(p 2 ), obtained by 
solving the quark DSE using the parameters in Eq. (|4l|), and in the chiral 
limit. It complements Fig. ^ by highlighting the qualitative difference between 
the behaviour of M(p 2 ) in the chiral limit and in the presence of explicit chiral 
symmetry breaking. In the latter case 



M{p 2 











p2 


r 


A 2 

. J1 QCD _ 





{ 1 + two loop} , 7„ 



12 



33 - 27V 



(42) 
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Figure 11: Dressed-quark mass function obtained in solving the quark DSE using the pa- 
rameters in Eq. (|tl|): u/d-quark (solid line); s-quark (long-dashed line); chiral limit (dashed 
line). The renormalisation point is £ m 20GeV. The intersection of the line M 2 (p 2 ) = p 2 
(dotted line) with each curve defines M E , the Euclidean constituent-quark mass. 



However, in the chiral limit the ultraviolet behaviour is 
largo-p 2 2?r 2 7 m (- (qq)°) 



M(p 2 ) 



In 



A 5 

QCD 



(43) 



where (qq}° is the renormalisation-point-independent vacuum quark conden- 
sate. Analysing the chiral limit solution yields 



(qq)° = (0.227 GeV) J , 



(44) 



which is a reliable means of determining (qq)° because corrections to Eq. (43) 
are suppressed by powers of Aq CD /£ 2 . 

Equation ( |36| ) defines the renormalisation-point-dependent vacuum quark 
condensate 



,_ , |C=19GeV 

w)c\ 



lim Z A {C,K)N C [ tr Dirac [5 A=0 (g)] 

A^oo J 



C=19GoV 



(45) 



and the calculatedEa value is (0.275 GeV) 3 . In the model it is straightforward 
to establish explicitly that mf- (qq)® = constant, independent of (, and hence 



24 



(46) 



which defines the renormalisation-point-independent current-quark masses un- 
ambiguously. From this and Eqs. (fid]), (|44|) and one obtains 

m u/d = 6.60 MeV , m s = 147 MeV . (47) 

Using the one-loop evolution in Eq. (|42] ) these values yield m^ d = 3.2 MeV and 
m'l = 72 MeV, which are within ~ 10% of the actual values in Eq. (41). This 



indicates that higher-loop corrections to the one-loop formulae provide contri- 
butions of < 10% at p 2 = ( 2 . The magnitude of the higher- loop contributions 
decreases with increasing p 2 . 

The renormalisation-point-invariant product in Eq. (|4^) also yields 

- (M>cl C =lGeV : = ( ln [VAqcd]) 7 "' (qq)° = (0-241 GeV) 3 , (48) 

which can be compared directly with the value ofj-the quark condensate em- 
ployed in contemporary phenomenological studies£U (0.236 ± 0.008 GeV) 3 . It 
is now straightforward to determine the accuracy of Eqs. ( |37| ) and (|8|). Using 
experimental values on the left-hand-side, one finds: 

(92.4 x 138. 5) 2 = (113 MeV) 4 cf. (Ill) 4 = 2 x 5.5 x 241 3 , (49) 
(113 x 495) 2 = (237 MeV) 4 cf. (206) 4 = (5.5 + 130) x 241 3 . (50) 

Hence, while it is good for the 7r, using the vacuum quark condensate for the 
kaon leads one to overestimate m^ GeV by 70%! 

4-1 Bethe-Salpeter Equation 

The explicit form of the Bethe-Salpeter equation consistent with Eq. ( p9] ) is 
= T H (k;P) + (51) 

-0^ G((k - q) 2 ) D^(k - q)^ S (q + )T H (q; P)S{q.)^ 7v = . 

With the Ansatz for Q{k 2 ) and the solution of the quark DSE for S(p), the 
kernel of the BSE is specified completely. It is a linear integral equation and 
solving to obtain Tff(k;P) and the bound state mass is a straightforward 
numerical problem. Having D^ u (k), S(p) and T^f(A;;P), the calculation of 
observables such as: the leptonic decay constant; meson charge radius, (qq) H ; 
and electromagnetic form factor, Fh(Q 2 ); etc., is possible. 
Using this exemplar one obtains, in MeV, 



f f° m f m VP = 1 ^ 2 fl p = 1 / 2 m r l p=U f 

Jtt J In ss j ss "Ik Ik " l k Jk (52) 



138.5 92.4 89.8 685 129 497 109 497 109 
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k z (GeV z ) 

Figure 12: Asymptotic behaviour of the Oth Chebyshev moments of the functions in the 
7r-meson Bethe-Salpeter amplitude: fTr°E 7r (k 2 ) (GeV, solid line); f n °i<V(A; 2 ) (dimension- 
less, long-dashed line); k 2 /„■ °G^-(fc 2 ) (dimensionless, dashed line); and k 2 °H 7v (k 2 ) (GeV, 
dot-dashed line). The fc 2 -dependence is identical to that of the chiral-limit quark mass func- 
tion, M(p 2 ), Eq. (pa). For other pseudoscalar mesons the momentum dependence of these 
functions is qualitatively the same, although the normalising magnitude differs. 



/° is the leptonic decay constant of a pseudoscalar meson in the chiral limit and 
"ss" denotes a fictitious pseudoscalar bound state. As emphasised by Eq. (|52|), 
observables are independent of the momentum partitioning parameter, r/p, 
when all the amplitudes in Eq. ( (30|) are retained in solving the BSE. 

One aspect of the Bethe-Salpeter amplitudes, which is particularly impor- 
tant in the calculation of electromagnetic form factors, is their behaviour at 
large-fc 2 . It is depicted in Fig. H2, where the zeroth Chebyshev moments are 



E H {k 2 ):=- f dp sin 2 (3 U {cos f3) E H (k 2 , k -P;P 2 ), 
n Jo 



(53) 



k ■ P := cos /3V k 2 P 2 , and similarly for the other functions. The momentum 
dependence of ^^(fc 2 ) at large-fc 2 is identical to that of the chiral-limit quark 
mass function, M(p 2 ) in Eq. ( |43| ) and characterises the form of the quark-quark 
interaction in the ultraviolet. Very importantly, the same is true of °F 7r (fc 2 ), 
fc 20 G w (fc 2 ) and fc 20 i/ 7r (fc 2 ), with analogous results for other mesons. This is 
critical, e.g., because it entailsn 2 ] that 

F^(q ) cx — 2 -f — ; (54) 



20 




i.e., q 2 F 7r (q 2 ) sa const., up to-calculable lng 2 -corrections. If one erroneously 
neglects F and G in IV, thenCJ q 4 F 7r (q 2 ) « const. 

5 Heavy Quarks 

The CKM matrix characterises the difference between the mass eigenstates (q) 
and weak eigenstates (</) in the Standard Model: 

*ud *us Vub 

V cd V cs V cb I I s | , (55) 
Vu Vtb 

and the Standard Model requires that the CKM matrix be unitary. The ma- 
trix elements V q Q are measurable in the semileptonic decay of a pseudoscalar 
meson: 

A(P Hl -> P H M = ^ V^q - 75)^ M^" 2 (pi.pj) , (56) 
where Gf is the Fermi weak-decay constant and the hadronic current is 

M^ Ph *( PuP2 ) := {PhMWiMP^M (57) 
= f+(t)( P i+ P2 )» + f-(t)q„, (58) 

with t := —q 2 - Hence accurate measurements and calculations of these decays 
can decide whether or not the Standard Model is complete. 

All the information about strong interaction effects in these processes is 
contained in the form factors, f±(t), appearing in the hadronic current, and 
their accurate estimation is essential to the extraction of V q Q from a measure- 
ment of a semileptonic decay rate: 

T(P Hl -» P H Jv) = (59) 



G 2 I 

mo 3 \V<lQ\ ~ 



— \ v «q\ 2 — / dt\f + (t)\ 2 [(t + -t)(t_ -t)f 2 



with t± := {rriH\ i 111 h 2 ) 2 an d neglecting the lepton mass. As with all form 
factors, the calculation of f± (t) requires a knowledge of the propagation char- 
acteristics of the mesons' constituents and the structure of the bound state. 
The DSEs are therefore ideally suited to this analysis. 



In Sec. 2.2 I introduced Cf, which measures whether a quark is light or 
heavy. If a quark is heavy then many simplifications ensue in the calculation 
of hadronic observables, and they can all be traced to being able to ignore the 
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momentum dependence of that quark's self energy. The 6-quark is certainly 
heavy but one cannot be as certain about the c-quark because its mass function 
exhibits significantly more momentum-dependence. Whether the c-quark is 
heavy is a qualitatively important quantitative question. It can be addressed 
by exploring the consequences of assuming 

S„(p) = - — y , S c (p) = - (60) 

where Mf ~ Mf. I will call this the "heavy-quark" limit. 

In this limit, when the other constituent is light, it is natural to work with 
i]p = 1 in Eq. ( |2^ ) since the total momentum of the bound state is dominated 
by the momentum of the heavy-quark. Introducing the heavy-meson velocity 
(v^) and binding energy (E) via 

P:=m H ^, v 2 = -l; M H :=M Q +E, (61) 
with Mh the heavy- meson mass, one findsEl at leading order in 1/M~h 

S( q + P)= 1 - 1 -^, (62) 
2 q ■ v — hi 

T H (q;P) = y/M^f H (q;P), (63) 

where the canonical normalisation condition for Th^F) is independent of 
M H . Using Eqs. @ and @ in Eq. © one obtains 

/^cxI/VMh. (64) 

This is obviously not correct for light quarks since < fx, and the experi- 
mental situation is depicted in Fig. ff3l It is clear from the figure that for the 
c-quark at least the 1/MQ-corrections to Eq. (Q) are significant. 

The mass formula of Sec. ^ has another corollary in the heavy-quark limit: 
using Eq. @ with Eqs. @ and @ yields 



Mh (x iriQ . (65) 



A model studyii shows this to be valid for toq > m s , which is confirmed 
by data, Fig. |lj. Although unsurprising when one considers that O ( m^)- 
corrections to Eq. (^) are large, as noted after Eq. (^p|), I observe that Eq. (|64| ) 
is not valid until at least toq > rh c . Equation (|6^) is then only valid because 
of cancellations between the 1/MQ-corrections to fn and (qq)f ■ 
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Figure 13: fjj as a function of the mass of the heaviest constituent, m q (I have included 
an extra, multiplicative factor of y/2, so that f n = 131 MeV, which is conventional in this 
application.) The solid line is ///(m 9 ) = const. /yjrh q , fitted to fg, while the dashed line 
assumes a 10% shift in fg from l/m q -corrections. 



5.1 Semileptonic Heavy — > Heavy Decays 

In impulse approximation the hadronic current describing Bf —> Df decays is 



r P Hl Ph 2 fm m s _ N c 
16tt 4 

tr [t H2 (k;-p 2 )Sq(k + p 2 )«V^ Q (fc + p 2 , k + pi)S Q {k + pi)T ' Hl (k;pi)S g > (fc)] , 



AC 1 a (Pi»Pa) = ^ / d ' k ( 66 ) 



where f H 2 (fc; —p-if '■= C^Th 2 (— fc; —p 2 )C, C — 7274, M* is the matrix trans- 
pose of M; and Vf^{kx,k 2 ) is the vector part of the dressed-quark-W-boson 
vertex. Using Eqs. ( |6^) and (|63|), and the approximation V^(ki,k 2 ) = 7^, 
which is valid in the heavy-quark limit, 

f± (*) = - - (w) ; (67) 

2 ^/m Df m Bf 

i.e., the form factors are determined by a single, universal function: £/(u>)@ 
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Fignre 14: Pseudoscalar meson mass as a function of the mass of the heaviest constituent, 
ra g l Only the it does not lie on the same straight line. Since rh s < M x , the mass-scale 
associated with DCSB, it is not too surprising that Eq. (^5| is valid for rhq > rh a : that is 
the domain on which explicit chiral symmetry breaking overwhelms DCSB. 



Employing an Ansatz [] for the heavy- meson Bethe-Salpeter amplitude 

f HlJ (k;pi) = 75 (l + \%1 ■ v) ±- ^(fc 2 ) , (68) 

where Kf is the fixed, canonical Bethe-Salpeter normalisation constant for 
t*H lf and <P(k 2 ) is any simple function that represents the heavy-meson as a 
finite-size, composite object; e.g., 

cj)(k 2 ) = cxp{-k 2 /A 2 ) , (69) 

with A a free fitting parameter, one obtains 

N f 1 1 f°° 

e Ladder-like truncations of the BSE are inadequate when the mass of one or both con- 
stituents becomes large; e.g., this truncation does not yield the Dirac equation when one 
constituent becomes infinitely massive. Pending improved BSE studies, an Ansatz is used. 



4(Z W ) 



<tUz w ) 



, (70) 
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f 

with o gjy the functions describing the propagation of the light-quark con- 



stituent, W = 1 + 2t(1 - t)(io - 1), z w = u - 2E v / u/W and f 

• v D 



w 



TOr + TO n — t 



(71) 



2to B/ to D/ 

The normalisation of the Bethe-Salpeter amplitude ensures automatically that 

0> = 1) = 1- (72) 



Experimentally measured violations of Eqs. (|67j) and (72) gauge the fidelity of 
the heavy-quark limit for physical processes. 



5.2 Semileptonic Heavy — > Light Decays 

The heavy-quark limit is not really helpful in studying the decays B — ► n, 
D — > K and D — ► 7r because there are only light-quarks in the final state and 
hence no useful expansion parameter. Therefore a theoretical description of 
these decays relies heavily on a good understanding of light quark propagation 
characteristics and the internal structure of light mesons. In this case, using 
Eqs. (K32) and (d3|), the impulse approximation to the hadronic current yieldsEj 



HiH 2 



V2 N c 



q fH 2 32** 



F q/ (t;E,m Hll m H2 ) 



where 



F q i(t]E,m Hll mH 2 ) = 

4 r 1 d 1 r 1 

7T y_i v/i -7 2 V 



(73) 



(74) 



dv 



u 2 du f(zi) £{z\) Wq' (7, v, u) 



with £ the pseudoscalar part of the light-meson Bcthc-Salpeter amplitude and 



W q ,(j,v,u) = 2i 



+ 



(75) 



1 - 



a u s (z 1 )a« v (z 2 ) + 



^The minimum physical value of w is ui m ; n = 1, which corresponds to maximum mo- 
mentum transfer with the final state meson at rest; the maximum value is w) ma x — 
(m| + m? D )/(2mg mjj ) = 1.6, which corresponds to maximum recoil of the final state 

meson with the charged lepton at rest. 
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°s(- 2; i) S (#2) +u^CTy(zi)cr| (z 2 ) + 

m Hl 

( Zl + uvM Hl )al{zi) a((z 2 ) - 2m 2 H2 r 2 ^-a( (z 2 ) 

where a'(x) := -£;0-(x) and 

zi = u 2 -2uvE, (76) 
Z2 = u 2 — 2u v [E — X) — m 2 H2 + 2i mH 2 7 u \/\ — v 2 , (77) 
X = (m Hl /2)[l + (m 2 H2 -t)/m 2 Hi ], (78) 

r = uv 7 !-^ 2 Vl- 7 2 - (79) 

Under the assumption of isospin symmetry a u also represents a <i-quark and, 
to illustrate Eq. (73|), the B° — ► 7r _ ^ + ^ decay is described by 

/,r 327H 

5.5 Calculated Transition Form Factors 

With the dressed-light-quark propagators and light-meson Bethe-Salpeter am- 
plitudes fixed completely in studiesu of ir- and if-meson properties the cal- 
culation of the heavy-meson transition form factors is straightforward. There 
are two free parameters: the binding energy E, introduced in Eq. (|6l]) and 
the width, A, of the heaw meson Bethe-Salpeter amplitude, introduced in 
Eq. (|69|). They were fixedtS at 

£ = 0.44GeV, A = 1.4 GeV = l/[0.14fm] (81) 

by requiring a best, weighted least-squares fiUo the three available lattice data 
pointsE£l for /?* and the experimental valueEZI for the B° — > n~£ + is branching 
ratio. In doing that the study used mg = 5.27 GeV and was constrained to 
yield fs = 0.17 GeV —which is the central value favoured in a recent analysis 
of lattice simulationscJ This procedure assumes only that the 6-quark is in the 
heavy-quark domain and yields Jb 3 = 0.18 GeV. 

The calculated form of /? 7r (t) is presented in Fig. [l5| A good interpolation 
of the result is provided by 

/+*(*) = 1 °w 58 2 . m ™n = 5.67 GeV . (82) 

1 - V^mon 
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0.0 10.0 20.0 

t ( GeV 2 ) 

Figure 15: Calculated form of f+ n (t). The solid line was obtained!^ assuming only that the 
fe-quark is heavy, the dashed line assumed the same of the c-quark. The data are the results 
of a lattice simulationESI and the light, short-dashed line is a vector dominance, monopole 
model: f+(t) = 0.46/(1 - t/m 2 B ,), m B , = 5.325 GeV. The light, dotted line is the phase 

space factor \f^{0)\ 2 [(t+ - t)(t_ - t)] 3/2 /(7rm s ) 3 in Eq. (||), which illustrates that the 
B — » irev branching ratio is determined primarily by the small-t .behaviour f^ n (t). The 
calculated branching ratio is 2.0 X 10~ 4 , cf. the experimental valuccll [1.8 ± 0.4] X 10 — 4 . 



This ptalue of m mon can be compared with that obtained in a fit to lattice 
data:ES m mon = 5.6 ± 0.3. The calculated value of ff n (0) = 0.46 is compared 
with its value obtained using a range of other theoretical tools in Table [j]. The 
uncertainty in f+*(0) is the major source of the large uncertainty in the CKM 
matrix element: V bu = 0.0032 ± 0.0012. 

The calculated form of £,(w), which characterises the semileptonic heavy- 
— > heavy-meson decays, is depicted in Fig. |l6|. It yields a value for p 2 := 
—$,'(w = 1) = 0.87 — 0.92, close to that obtained with a linear fitting formal 

t fit. The 
however, 



however, £,(w) has significant curvature and deviates quickly from t 
curvature is, in fact, very well matched to that of the nonlinear fit 
the value of p 2 reported in that case is very different from the calculated value. 
The derivation of the formula for £,(w) assumes that the heavy-quark limit is 
valid not only for the &-quark but also for the c-quark. These results therefore 
suggest that the latter assumption is only accurate to approximately 20%. 

The calculated form of f+ K (t) is depicted in Fig. [Tt[. The i-dependence is 
well-approximated by a monopole fit. The calculated pvalue of f+ K (0) = 0.62 
is approximately 15% less than the experimental valuetl That is also a gauge of 
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Reference 


/f-(0) 


Ref. f§ 


0.46 


Dispersion rel 
Quark Model 
Quark Model 
Quark Model 

Light-Cone Si 

Quark Confin 
Quark Confin 


ai 

30 
31 
32 

IE 

er 
er 


ions 1^1 

i Rules S 

nent Model 2f 
nent ModelSH 


0.18 -> 0.49 
0.33 ±0.06 
0.21 ±0.02 
0.29 

f 0.29 direct 

| 0.44 pole dominance 

0.6 

0.53 



Table 1: Comparison of the calculated result for f^(0) with a values obtained usmo-otiifir 
theoretical tools. More extensive and complementary lists are presented elsewhereEj'E3c3 



the size of l/M c -corrections, which are expected to reduce the value of the D- 
and _D s -meson leptonic decay constants calculated in the heavy-quark limit: 
f D = 285 MeV, f Ds = 298 MeV. A 15% reduction yields f D = 0.,24 GeV and 
So, = 0.26 GeV 
with experiment 



lues which are consistent with lattice estimatespj the latter 
and both with the scaling law illustrated in Fig. [l^. 
On the kinematically accessible domain, < t < (mo — m^) 2 , the calcu- 
lated form of S+ n {t) is well described by the monopole fit 

0.716 



/+*(*) = 



1 — t/m 



i 

mon 



= 2.15 GeV. 



(83) 



A naive vector .meson dominance assumption would employ m D 
2.0 GeV. UsingE! \V cd /V cs \ 2 = 0.051 ± 0.0020 



rriD* 



Rtt 



Br(D -> -kZv) 
Br(D -> Klv) 



2A7 



V cd 



V r , 



0.13. 



(84) 



The bulk of the l/M c -corrections cancel in this ratio. Experimental!; 



y§§§ 



Rji 
Rtt = 2 



Br(D° 



it e 



Br(D° 
Br{D+ 



= 0.11™ ±0.1 



= 0.17 ±0.05 ±0.03 



(85) 



(86) 



~Br(D+ -> K a e+v e ) 

If one assumes single-pole, D* and D* vector meson dominance for the t- 
dependence of the form factors S+ and S+ > one obtains the simple formula 



Rtt = 1.97 



S?*(0) 



S° K (0) 



Vcd 



(87) 
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Figure lj&ri Calculated form of £(u>) cf. recent experimental analyses. The solid line was 
obtainecO assuming only that the 6-quark is_heavy, the dash-dot line assumed the sa mp o f 
the c-quark. Experiment: data points - Ref. |pj| ; short-dashed line - linear fit from Ref. |2S i 
f(iu) = 1-p 2 (w-1), p 2 = 0.91 ± 0.15 ±0.16 ; long-dashed line - nonlinear fit from Ref. 38 , 
£{w) = [2/(w + 1)] exp [(1 - 2p 2 ) (w - l)/(w + 1)] , p 2 = 1.53 ± 0.36 ± 0.14. The two 
light, dotted lines are this nonlinear fit evaluated with the extreme values of p 2 : upper line, 
p 2 = 1.17 and lower line, p 2 = 1.89. 



That approach was usee 
1.3 ±0.2 ±0.1 from Eqs. 



= 1.0 



J to estimate f£*(0)/f? K (0) ^ 
(H) and @. This calculation^ yields" 



\ll ± 0.04 or 



/f K (0) 



1.16. 



It must be noted that Ref. |19f| explicitly did not assume vector meson dom- 
inance. The calculated results reflect only the importance and influence of the 
dressed-quark and -gluon substructure of the heavy mesons. That substruc- 
ture is manifest in the dressed propagators and bound state amplitudes, which 
fully determine the value of every calculated quantity. Explicit vector meson 
contributions would appear as pole terms in V^ 1 ^ 2 (fci, fe), which are excluded 
in the Ansatz: VH (k\,k2) — 7^- That simple-pole Ansatze provide efficacious 
interpolations of the calculated results on the accessible kinematic domain is 
not surprising, given that the form factor must rise slowly away from its value 
at t = and the heavy meson mass provides a dominant intrinsic scale, which 
is only modified slightly by the scale in the light-quark propagators and meson 
bound state amplitudes. Similar observations are true in the calculation of the 
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Figure 17: Calculated form of f+ K (q 2 ): the solid line was obtaine assuming only that the 
6-quark is heavy, the dashed line assumed the same of the c-quark. The light, short-dashed 
line is a vector dominance, monopole model: f+(q 2 ) = 0.74/(1 — q 2 /m 2 „ ), m D * = 2.11GeV. 



pion form factorE3e3 

In this section I illustrated a heavy-quark limit of the DSEs, based on 
the observation that the mass function of heavy-quarks evolves slowly with 
momentum, and the manner in which it can be used to describe the properties 
of heavy-mesons. They are much like light-mesons: bound states of finite 
extent, with dressed-quark constituents. The analysis of B — > D, D — > K and 
D — > 7r transitions indicates that the heavy-quark limit is accurate to within 
15-20% for the c-quark. A significant feature is the correlation of heavy — > 
heavy and heavy — > light transitions and their correlation with light meson 
observables, which are dominated by effects such as confinement and DCSB. 



6 Finite T and fi 

Nonperturbative methods are necessary to study the transition to a QGP, 
which is characterised by qualitative changes in order parameters such as the 
quark condensate. One widely used approach is the numerical simulation of 
finite-T lattice-QCD, wi 
tensive efforts since then 



the first simulations in the early eighties and ex- 
They have contributed considerably to the current 



3G 



understanding of the nature of the QGP. The commonly used quenched approx- 
imation is inadequate for studying the phase diagram of finite-T QCD because 
the details of the transition depend sensitively on the number of active (light) 
flavours. It is therefore necessary to include the fermion determinant. 

That is even more important in the presence of p,, which modifies the 
fermion piece of the Euclidean action: 7 • d + m — > 7 • d — 74/i + rn, introducing 
an explicit imaginary part to the fermion determinant. With the p ^ lattice- 
QCD action being complex, the study of finite density is significantly more 
difficult than that of finite temperature. Simulations that ignore the fewnion 
determinant encounter a forbidden region, which begins at p = m,r/2p3 and 
since — ► in the chiral limit this is a serious limitation, preventing a reliable 
study of chiral symmetry restoration. i-The phase of the fermion determinant 
is essential in eliminating this artefact ej 

The contemporary application of DSEs at finite temperature rand chemical 
potential is a straightforward extension of the T = = fi studiesc3 The direct 
approach is to develop a finite-T extension of Ansdtze for the dressed-gluon 
propagator. The quark DSE can then be solved and, having the dressed- 
quark and -gluon propagators, the response of bound states to increases in 
T and /x can be calculated. As a nonperturbative approach that allows the 
simultaneous study of confinement and DCSB, the DSEs have a significant 
overlap with lattice simulations: each quantity that can be estimated using 
lattice simulations can also be calculated using DSEs. That means they can be 
used to check the lattice simulations, and importantly, that lattice simulations 
can be used to constrain their model-dependent aspects. Once agreement is 
obtained on the common domain, the DSEs can be used to explore phenomena 
presently inaccessible to lattice simulations. 

6.1 Finite-(T, fi) Quark DSE 

The renormalised dressed-quark propagator at finite- (T, fi) has the form 



where Cu k := tu k + ip with u> k = (2fc + 1) 7rT, fee Z, the fermion Matsubara 
frequencies. The complex scalar functions: A(p, B(p,ui k ) and C{p,u)k) 
satisfy: T{p, cuk)* = Fifi, &-k-i) i ? = ^,fi,C, and although not explife 
itly indicated they are functions only of |p| 2 and w\. Using the Matsubarau3 
(or imaginary-time) formalism 0(4) covariance, the Euclidean realisation of 



S(p,u>k) 



1 



(89) 



17 • pA(p 7 u; k ) + ij 4 u; k C(p,u} k ) + B(p,u) k ) 
~ij ■ p a a (p, oj k ) - ij4 Lu k <j c (p, u k ) +<Jb (p, 0Jk) 



(90) 
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Poincare covariance, is broken to 0(3) and only systems in equilibrium can be 
studied. 

S(p,£jk) satisfies the DSE 

S~ x {p, u k ) = Z£ ij ' P + z 2 («74 Wfe + m bm ) + uj k ) , (91) 

where the regularised self energy is 

E'(p,u) k ) = ■ pY,' A (p, Q k ) +«74Wfc£' c (p,u)fc) +Y,' B (p,u) k ) , (92) 

Jl.q 6 4 

(93) 

It ■■= T ZZ-oo I A ^d V A := -(Zf/p 2 ^J ■ P, V B := Z u V c := 

— (Zi/Qk)ij4- The finite- (T, fi), Landau-gauge dressed-gluon propagator in 
the kernel has the form 

g 2 D^(p, n) = p^(p, n) A F (p, n) + p?m a g ( p , n) , (94) 

{0; /i and/or v = 4, 

^ = ^ = 1,2,3 - ( 95 ) 

with Pj y (p) + P^{p,p 4 ) = -p M P./(El=i PqPo); M,^ = 1,...,4. 

Whether or not (T, //) are nonzero, in studying confinement one cannot 
assume that the analytic structure of a dressed propagator is the same as 
that of the free particle propagator: it must be determined dynamically. The 
pk '■= (p, (I>fc)-dependence of A and C is qualitatively important since it can 
conspire with that of B to eliminate free-particle poles in the dressed-quark 
propagatorCJ In that case the propagator does not have a Lehmann represen- 
tation so that, in general, the Matsubara sum cannot be evaluated analytically. 
More importantly, it further complicates a real-time formulation of the finite 
temperature theory making the study of nonequilibrium thermodynamics a 
very challenging problem. In addition, the p^-dependence of A and C can be a 
crucial factor in determining the behaviour of bulk thermodynamic quantities 
such as the pressure and entropy, being responsible for these quantities reach- 
ing their respective ultrarelativistic limits only for very large values of T and 
/.t. It is therefore important in any DSE study to retain A(pk) and C(pk), and 
their dependence on pk . 
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6.2 Phase Transitions and Order Parameters 

One order parameter for the chiral symmetry restoration transition is the quark 
condensate, defined via the renormalised dressed-quark propagator, Eq. (|36|). 
An equivalent order parameter is 

X := Re£? (p = 0,u) ) , (96) 

which makes it clear that the zeroth Matsubara mode determines the charac- 
ter of the chiral phase transition. An order parameter for confinement, valid 
for both light- and heavy-quarks, was introduced in Ref. Q. It is a single, 
quantitative measure of whether ox. not a Schwinger function has a Lehmann 
representation, and has been usedEJ to striking effect in QED3. 



6.3 Study at (T ^0,^ = 0) 

Deconfinement and chiral symmetry restoration at T 7^ have been studiedEl 
in a DSE-model of two-light-flavour QCD. The quark DSE was solved using 
a one-parameter model dressed-gluon propagator, which, provides a good de- 
scription of 7r and p-meson observables at T = = /.ilj The transitions are 
coincident and second-order at a critical temperature of T c w 150 MeV, with 
the same estimated critical exponents: j3 — 0.33 ± 0.03. Both and /„• 
are insensitive to T for T < 0.7 T c . The evolution with T is so slow that 
even at T = 0.9 T c there is only a 20% suppression of T^^gy However, for T 
very near to T c the pion mass increases substantially, as thermal fluctuations 
overwhelm quark-antiquark attraction in the pseudoscalar channel, until, at 
T = T c , — ► and there is no bound state. These results paifirmQ those of 
contemporary numerical simulations of finite-T lattice-QCDu3 



6.4 Complementary study at (T = 0, \i ^ 0) 



The behaviour of this model at-it 7^ has also been explored. I 
simple dressed-gluon propagatortiil 



Using the 



G(k 2 



16 2 

7T 

9 



4ir 2 m 2 t 6 4 {k) 



-[fe 2 /(4m?)] 



k 2 



(97) 

9 The lattice and DSE cstimat*s-ipi-i(3,Y , however, do not survive more exhaustive studies 
and the most recent analyse faufl suggest that in DSE models whose long-range part is 
described by a (5-function singularity the chiral symmetry restoration transition at finitc-T 
is described by a mean-field value of /3. 
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where m f = 0.69 GeV= l/[0.29fm] is a mass-scale that marks the boundary 
between the perturbative and nonperturbative domains, the quark DSE was 
solved in rainbow approximation. The solution has the form 

S(P[»]) ■= -ii'PVAiPw) -*74<«V]<rc(P|>]) +o- b (pm), (98) 
where p [fl] := (p, aj [fl] ), with uj [fl] := p A + ifi 



As elucidated in Sec. 6J5, there are two distinct types of solution: a Nambu- 
Goldstone mode with confinement and DCSB characterised by <jb ^ 0; and a 
deconfined, chirally-symmetric Wigner-Weyl mode characterised by ctb = 0. 
The possibility of a phase transition between the two modes is explored by 
calculating the relative stability of the different phases, which is measured by 
the difference between their tree-level auxiliary-field effective-actions: 

W : = / 7^1 (99) 



2Af / Af c vrv J (2tt 



In 



IpTAg+o^Cg + Bg 



+ \P\ 2 (°A ~ &a ) + wf^j (a Co - CTCo) 



where A and C represent the solution of Eq. ( |9l"| ) obtained when Bq = 0; 
i.e., when DCSB is absent— This solution exists for all ix. B(fx) defines a fi- 
dependent "bag constant" It is positive when the Nambu-Goldstone phase 
is dynamically favoured; i.e., has the highest pressure, and becomes negative 
when the Wigner pressure becomes larger. Hence the critical chemical potential 
is the zero of B(fi), which is fi c = 375 MeV. The abrupt switch from the Nambu- 
Goldstone to the Wigner-Weyl mode signals a first order transition. 

The chiral order parameter increases with increasing chemical potential 
up to jU c , with X(/i c )/X(0) « 1.2, whereas the deconfincment order param- 
eter, k(ji), is insensitive to increasing /i. At /i c they both drop immediately 
and discontinuously to zero, as expected of coincident, first-order phase tran- 
sitions. The increase of X with /i is a necessary consequence of the momentum 
dependence of the scalar piece of the quark self energy, Z?(p^])L3 The vacuum 
quark condensate behaves in qualitatively the same manner as X . 

Even though the chiral order parameter increases with /i, decreases 
slowly as /i increases, with m 7r (/i ~ 0.7/i c )/m 7r (0) ~ 0.94. At this point 
m v begins to increase although, for /i < /i c , m^/i) does not exceed m w (Q), 
which precludes pion condensation. The behaviour of m w results from mutually 
compensating increases in f% and (m^qq) f) „ . f n is insensitive to the chemical 
potential until \i rs 0.7 (i c , when it increases sharply so that f v (/i~ (/i = 
0) ~ 1.25. The relative insensitivity of and / w to changes in /i, until very 
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near p c , mirrors the behaviour of these observables at finite-TEj For example, 
it leads only to a 14% change in T^—^^ at fi w 0.9 p c : an increase in this 
case. The universal scaling conjecture of Ref. |58) is inconsistent with the 
anticorrelation observed between the /z-dependence of f v and m ff . 

6.4-1- /i-dependence anticorrelated with T-dependence. Comparing the 
/i-dependence of f w and m n with their T-dependence, one observes an anticor- 
relation; e.g., at fi = 0, f„ falls continuously to zero as T is increased towards 
T c fts 150 MeV. This is a necessary consequence of the momentum-dependence 
of the quark self energy. In calculating these observables one obtains expres- 
sions for m 2 or / 2 and thus the natural dimension is mass-squared. Therefore 
their behaviour at finite T and \x is determined by 

Re(a;f M] )~[7r 2 r 2 - M 2 ], (100) 

where the T-dependence arises from the introduction of the fermion Matsubara 
frequency: p4 — > (2k + l)irT. Hence_syhen such a quantity decreases with T it 
will increase with p, and vice-versaE3 

6-4-2: (—(qq)) increases with p. The confined-quark vacuum consists of 
quark-antiquark pairs correlated in a scalar condensate and increasing ji in- 
creases the scalar density: (— (qq) ) . This is an expected consequence of confine- 
ment, which entails that each additional quark must be locally paired with an 
antiquark thereby increasing the density of condensate pairs as p is increased. 
For this reason, as long as fj, < fi c , there is no excess of particles over an= 
tiparticles in the vacuum and hence the baryon number density remains zeropO 
i.e., p U Q d = , (j, < fjb c . This is just the statement that quark-antiquark pairs 
confined in the condensate do not contribute to the baryon number density. 

6-4-3: The core_of neutron stars. The vacuum quark pressure, P u+d [fj], 
can be calculatedE3 After deconfinement it increases rapidly, as the condensate 
"breaks-up" , and an excess of quarks over antiquarks develops. The baryon- 
number density, p u ^ d — (l/"i)dP u+d / 'dfi, also increases rapidly, with 

pl+ d (li^2p c )~Z Po , (101) 

where po = 0.16 fm~ 3 is the equilibrium density of nuclear matter. For compax^ 
ison, the central core density expected in a 1.4 M@ neutron star is 3.6-4.1 poEi 
Finally, at p ~ 5/i c , the quark pressure saturates the ultrarelativistic limit: 
pu+d = y±/(2n 2 ) : and there is a simple relation between baryon-density and 
chemical-potential: 

P U B F+dF M = \Ki Va*>5 Mc , (102) 

6 IT 
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so that p U j[ F (5/i c ) ~ 350 po- Thus the quark pressure in the deconfincd 
domain overwhelms any finite, additive contribution of hadrons to the equation 
of state. That was anticipated in Ref. |55f| where the hadron contribution was 
neglected. This discussion suggests that a QGP is likely to exist in the core of 
dense neutron stars. 



6.5 Simultaneous study of (T 7^ 0, /i 7^ 0) 

This is a difficult problem and the most complete studyEll to date employs 
a simple Ansatz for the dressed-gluon propagator that exhibits the infrared 
enhancement suggested by Ref. Il7) : 



fD,M n*) = [s^ prf^j ^ y 6k0 ' (103) 

with rik = 2kirT, the boson Matsubara frequency. As an infrared-dominant 
model that does not represent well the behaviour of D MJ ,(p, fife) away from 
\p\ 2 + fi| w 0, some model-dependent artefacts arise. However, there is signif- 
icant merit in its simplicity and, since the artefacts are easily identified, the 
model remains useful as a means of elucidating many of the qualitative features 
of more sophisticated Ansatze. _ . 
With this model, using the rainbow approximation, the quark DSE ia!3 

S ,_1 (p,w fc ) = S^ipT&k) + ^JvSiPi^^jv . (104) 



A simplicity inherent in Eq. (103) is now apparent: it allows the reduction of 
an integral equation to an algebraic equation, in whose solution many of the 
qualitative features of more sophisticated models are manifest. 



In the chiral limit Eq. (104) reduces to a quadratic equation for B(pk), 
which has two qualitatively distinct solutions. The Nambu-Goldstone solution, 
with 



5fc)^^^' R «<T 
, otherwise 



(105) 

C{ ^ ) = \ i f, , ' ( 106 ) 



2 



Re(pl) < £ 



otherwise 



describes a phase of this model in which: 1) chiral symmetry is dynamically 
broken, because one has a nonzero quark mass-function, B(pk), in the absence 
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of a current-quark mass; and 2) the dressed-quarks are confined, because the 
propagator described by these functions does not have a Lehmann representa- 
tion. The alternative Wigner solution, for which 



B(p k ) = , C(p k ) = \ (l + ^1 + |0 , (107) 

describes a phase of the model with neither DCSB nor confinement. Here 
the relative stability of the different phases is measured by a (T, ^-dependent 
vacuum pressure difference; i.e., a (T, /^-dependent bag constant: B(T,p). 

B(T, p) = defines the phase boundary, and the deconfinement and chiral 
symmetry restoration transitions are coincident. For p = the transition is 
second order and the critical temperature is Tj? = 0. 159 77, which using the 
value of 77 = 1.06 GcV obtained by fitting the ir and p masses corresponds to 



T c ° = 0.170 GeV. This is only 12% larger than the value reported in Sec. 6.3 
and the order of the transition is the same. For any p ^ the transition is 
first-order. For T = the critical chemi cal p otential is p° c = 0.3 GeV, which 
is « 30% smaller than the result in Sec. |6.4| . The discontinuity in the order 
parameter vanishes at p = 0, and this defines a tricritical point. 

In this model the quark pressure, P q , is calculated easily. Confinement 
means that P q = in the confined domain. In the deconfined domain it 
approaches the ultrarelativistic, free particle limit, Ptjr, at large values of 
T and p. However, the approach to this limit is slow. For example, at 

simi- 
This 



T ~ O.377 ~ 2T c, or p ~ 7) ~ 3p° c , P q is only 0.5 P UR . A qualitative 
lar result is observed in numerical simulations of finitc-T latticc-QCD 
feature results from the persistence of momentum dependent modifications of 
the quark propagator into the deconfined domain, as evident with C ^ 1 in 



4E 



Eq. (107). It predicts a "mirroring" of finite-T behaviour in the /z-dependence 
of the bulk thermodynamic quantities, as illustrated in Fig. [l^. 



6.6 7r and p properties 

The model introduced in the last section has also been usec@ to study the 
(T, /^-dependence of meson properties, and to elucidate other features of mod- 
els that employ a more sophisticated Ansatz for the dressed-gluon propagator. 
For example, the vacuum quark condensate is given by the simple expression 

- (qq) = rf ^ r dy y 2 Re (J\-y 2 -u>?) , (108) 

T = T/r), p, — p/rj; l max is the largest value of I for which clf max < (1/4) + p 2 
and this also specifies w/ max , A 2 = uif — u>f, pi = (y, Qi+ifi). At T = = p. 



43 



Figure 18: The quark pressure, P q (T,/i), normalised to the free, massless (or Ultra- 
Relativistic) result. 



{-(m)) = V 3 /(80tt 2 ) = (0.11 ?7) 3 . Obvious from Eq. fllQ8| ) is that (-(qq)) 
decreases continuously to zero with T but increases with /i, up to a critical value 
of Hc(T) when it drops discontinuously to zero: behaviour in agreement with 
that reported in Sees. 6.3 and 5.4. The vacuum rearrangement emphasised in 



Sec. 6.4.2 is manifest in the behaviour of the necess arily - momentum-dependent 
scalar part of the quark self energy, B[pk) in Eq. (105). 

The leptonic decay constant is also given by a simple expression in the 
chiral limit: 



A? 



1=0 



1 + 4/T -4u>f 



A? 



(109) 



As anticipated, the co mbin ation /i 2 — uif, emphasised in Sec. 6.4.1 , cha racterises 
the behaviour of Eqs. (108) and (109). Without calculation, Eq. (109) indicates 
that will decrease with T and increase with /z, which provides a simple 
elucidation of the results in Sees. jfO] and [^4| 

The (T, /^-response of meson masses is determined by the ladder BSE 



2 

T M (pk;Pe) = -\ Re { 7 „ S(pi + \Pt) r M (pi! Pi) S(p t - \P t ) 7m 



) , (no) 
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1.0 



M pt ,T = 

M„, T = 

M pt ,T = 0.1r| 



a 0.6 



B 0.4 



0.2 - 



0.0 



0.0 



0.1 



0.2 



0.3 



Figure 19: M p + = M u + and as a function of /i for T = 0, 0.1. On the scale of this 
figure, M„ is insensitive to this variation of T. The current-quark mass is m = 0.011 r?, 
which for r) = 1.06 GeV yields M p+ = 770 MeV and = 140 MeV at T = = fi. 



where Pg := (P, Qe), with the bound state mass obtained by considering p = o- 
In this truncation the u>- and p-mesons are degenerate. 
The equation admits a solution for the 7r-meson with 



(111) 



and the calculated (T, /^-dependence of the mass is depicted in Fig. [L9[ 

For the p-meson the solution has two components: one longitudinal and 
one transverse to P. The solution of the BSE has the form 




\P\ 2 



(112) 



where 9 p+ labels the longitudinal and 9 P - the transverse solution. The eigen- 
value equation obtained from Eq. (110) for the bound state mass, M p ±, is 



if 



Re 



± CO, 



0+ 



cry (u, 



-K± 

4 P ± 



Th e ma ss of the transverse component is obtain ed w ith [— ujq + — (1/4)M^_] 
in Eq. (113). Using the chiral-limit solutions, Eq. (105), one obtains immedi- 
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ately that 

M p _ = — rj 2 , independent of T and p. (H4) 

Even for nonzero current-quark mass, M p _ changes by less than 1% as T and 
/x are increased from zero toward their critical values. Its insensitivity is con- 
sistent with the absence of a constant mass-shift in the transverse polarisation 
tensor for a gauge-boson. 

For the longitudinal component one obtains in the chiral limit: 

Af 2 + = \r? - 4( M 2 - . 2 T 2 ) . (115) 

The characteristic combination p 2 — 7r 2 T 2 again indicates the anticorrelation 
between the response of M p+ to T and its response to p, and, like a gauge- 
boson Debye mass, that M p+ rises linearly with T 2 for p = 0. The m ^ 
solution of Eq. Qll3| ) for the longitudinal component is plotted in Fig. [ll| M p+ 
increases with increasing T and decreases as p increases. 

Equation ( |113| ) can also be applied to the </>-meson. The transverse compo- 
nent is insensitive to T and p, and the behaviour of the longitudinal mass, M^ + , 
is qualitatively the same as that of the p-meson: it increases with T and de- 
creases with p. Using r\ = 1.06 GeV, this simple model yields M$± = 1.02 GeV 
for m s = 180 MeV at T = = p. 

In a 2-flavour, free-quark gas at T = the baryon number density is 
Pb = 2/i 3 /(37r 2 ) , by which gauge nuclear matter density, po — 0.16 fm -3 , 
corresponds to p = po := 260 MeV = 0.245^. At this chemical potential the 
algebraic model yields 

M p+ (/i ) « 0.75M p+ (p = 0) , M$ + (p ) w O.85M 0+ (/i = 0) . (116) 



The study summarised in Sec. 6.4 indicates that a better representation of the 
ultraviolet behaviour of the dressed-gluon propagator expands the horizontal 
scale in Fig. jl9L with the critical chemical potential increased by 25%. This 
suggests that a more realistic estimate is obtained by evaluating the mass at 
Pq = 0.20 r], which yields 

M p+ {p' ) « 0.85M p+ (p = 0) , M$ + (p' ) pb 0.90A^+(/i = 0) ; (117) 



a small, quantitative modification. The difference between Eqs. (116) and 



(117) is a measure of the theoretical uncertainty in the estimates in each case. 
Pursuing this suggestion further, p — \/2p' Q , corresponds to 2p , at which 
point M u+ = M p+ w 0.72 Mp+(jJL = 0) and Af 0+ w 0.85 M^+(/i = 0), while at 
the T — critical chemical potential, which corresponds to approximately Zpa 
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in Sec. K$, M u+ = M p+ » 0.65Afp+(/j = 0) and M 0+ » 0.80 M 0+ (/i = 0). 



These are the maximum possible reductions in the meson masses. 

This simple model preserves the momentum-dependence of gluon and quark 
dressing, which is an important qualitative feature of more sophisticated stud- 
ies. Its simplicity means that many of the consequences of that dressing can be 
demonstrated algebraically. For example, it elucidates the origin of an anticor- 
relation, highlighted in Sec. 6.4.1 and found for a range of quantities, between 
their response to increasing T and that to increasing fx. That makes clear why 
the transition to a QGP is second order with increasing T and first order with 
/i. Further, in providing an algebraic explanation of why the (T, /^-dependence 
of (—(qq)) and / w is opposite to that observed for M p+ , it emphasises that the 
scaling law conjectured in Ref. []58j is inconsistent with many actual calcula- 
tions that preserve the global symmetries of QCD. 



7 Some Final Remarks 

There are many challenges in hadronic physics and I have only presented an 
outline here. What should be clear, however, is that this field is prescribed by 
the need for nonperturbative methods and models. The models are necessary 
to propagate the qualitatively robust results of difficult nonperturbative studies 
into that realm where comparison with the current and future generation of 
experiments is possible. 

The DSEs serve both purposes: they provide a nonperturbative tool and a 
model framework, and have been used widely and efficaciously. In addition to 
the exemplars described herein, there are many other recent studies. Some are 
reviewed in Ref. [Q , others are axalculation of the cross sections for diffractrue, 
vector meson electroproductionjEI the electric dipole moment of the /?-mesonp3 
and an exploration of 77-77' mixingS 

There are two key phenomenological aspects of contemporary DSE appli- 
cations: modelling the infrared behaviour of the gluon propagator and trun- 
cating the kernel of the Bethe-Salpeter equation. They introduce a model- 
dependence, which is restricted to the infrared because the weak-coupling ex- 
pansion reproduces perturbation theory and hence the k 2 > 1 GeV 2 behaviour 
is fixed in both cases. 

As always, one should maintain a constructive scepticism about the fidelity 
of model input. Its impact must be gauged; e.g., by exploring and exploiting 
the constraints that Ward Identities and Slavnov- Taylor identities can nipvide 
in the theory. That approach has been particularly fruitful in QED^3 and 
already in the development]^ of a systematic truncation procedure for the 
kernel of the BSE in QCD. 
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Such checks are a useful adjunct. However, widespread phenomenological 
applications remain a key. Pushing a single framework, successful on a large 
domain, to its limits is a very effective means of identifying which elements of 
the approach need improvement, and an excellent way to focus resources on 
those sites requiring further, fundamental developments. 
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